astro-ph /0403597| 



Hyperbolic Universes with a Horned Topology and 
the CMB Anisotropy 

■ Ralf Aurich^, Sven Lustig^, Frank Steiner^ and Holger Then^ 



' ^Abteilung Theoretische Physik, Universitat Ulm, 

^ . Albert-Einstein- Allee 11, D-89069 Ulm, Germany 

^School of Mathematics, University of Bristol, 
University Walk, Bristol, BS8 ITW, United Kingdom 



O 



Abstract. We analyse the anisotropy of the cosmic microwave background (CMB) in 
hyperbolic universes possessing a non-trivial topology with a fundamental cell having 
' an infinitely long horn. The aim of this paper is twofold. On the one hand, we show 

, that the horned topology does not lead to a flat spot in the CMB sky maps in the 

' direction of the horn as stated in the literature. On the other, we demonstrate that 

l/^ I a horned topology having a finite volume could explain the suppression of the lower 

I multipoles in the CMB anisotropy as observed by COBE and WMAP. 

o 

^ ; PACS numbers: 98.80.-k, 98.70.Vc, 98.80.Es 
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1. Introduction 



A fundamental problem in cosmology is the large-scale geometry of the Universe, in 
particular its spatial curvature and topology. Since the Einstein gravitational field 
^ ' equations are differential equations, they determine the local properties of space-time, 
c3 ] but not the global structure of the Universe at large. In the so-called concordance 
model of cosmology it is assumed that the Universe is at large scales spatially fiat and 
possesses the trivial topology, implying that it has infinite volume. In the framework 
of inflationary scenarios, these properties are supposed to be determined by the initial 
conditions at the Big Bang. However, since we are lacking a theory of quantum gravity, 
the initial conditions cannot be derived from first principles. Instead, one can analyse 
the recent astronomical data, in particular on the cosmic microwave background (CMB) 
radiation, in order to deduce the curvature and the topology of the Universe at large 
scales. 

In 1992, COBE jjj made the spectacular discovery of the temperature fluctuations 
of the CMB and thus provided important clues about the early Universe and its time 
evolution. Expanding the observed temperature fluctuations ST{n) across the microwave 
sky into spherical harmonics Yj^{n), yields the expansion coefficients aim which in turn 
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lead to the multipole moments 
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and the angular power spectrum 6T^ := /(/+l)Q/(27r). In particular, COBE Ij detected 
in the angular power spectrum of the CMB a low quadrupole moment C2 corresponding 
to a strange suppression of power on large angular scales. It was soon realized that 
although the standard cosmological models are in agreement with the CMB anisotropy 
on small and medium scales, they fail to match the loss of power on large angular scales, 
especially those corresponding to the quadrupole moment. This observation was one of 
the motivations to study non-trivial topologies (see the reviews [21 IS])- In particular, 
compact hyperbolic universes were studied by several authors jH El El El El El UHI Ulj • For 
these compact hyperbolic universes, it was shown that the non-trivial topology leads 
indeed to a suppression of Ci for small values of /. Thus the observed loss of power 
on large angular scales was interpreted as a clear hint to a non-trivial topology of our 
Universe. 

The first findings of NASA's explorer mission "Wilkinson Microwave Anisotropy 
Probe" (WMAP) J2] have tremendously increased our knowledge of the temperature 
fluctuations of the CMB, since WMAP has measured the anisotropy of the CMB 
radiation over the full sky with high accuracy. The WMAP data confirm not only 
COBE's measurement of the low quadrupole moment C2, but display in the temperature 
(auto-) correlation function C{-^) {h ■ h' = cos'd) 



very weak correlations at wide angles, 70° ^ ^ ^ 150°, see the dashed curve in figure 
m (Note that the monopole and dipole are not included in the sum (j21).) At the largest 
angles, d > 160°, the WMAP-data display even a "correlation hole", i. e. negative values 
of CiJ}). In figure m we also show as a dotted curve the theoretical prediction according 
to the concordance model using the best-fit values for the cosmological parameters as 
obtained by WMAP ^21- It is seen that the concordance model, i. e. the best-fit ACDM 
model, does not reproduce the experimentally observed suppression at ^9 > 60° and the 
observed correlation hole. 

Since the correlation function C{d) emphasizes large angular scales and thus the 
low /-range, it is an ideal indicator function to search for a fingerprint of a possible 
non-trivial topology of the Universe. 

The evidence for a low quadrupole and missing power on large scales is, however, 
currently heavily discussed. The low value of the quadrupole obtained by the WMAP 
team ^2] has been criticized by the way how the foregrounds 112] are taken into account 
which are mainly caused by free-free, synchrotron and dust emission. Alternative 
reconstructions of the true cosmological signal from the five frequency bands measured 
by WMAP are discussed in jTH HHj which lead to a higher quadrupole moment. The 
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Figure 1. The correlation function C(i9) from the WMAP-data (dashed curve) in 
comparison with the concordance model (dotted curve) . 

observation that the plane of the quadrupole and two of the three planes of the octopoles 
are aligned towards the ecliptic might be seen as a hint of an unknown source or sink 
of CMB radiation in the outer solar system or as an unrecognized systematic 
The influence of the approximation of the likelihood function used for the angular 
power spectra analysis is investigated in where for the lowest multipoles an exact 
likelihood function estimation is presented. The statistical significance of the low value 
of the quadrupole moment is also discussed in [TH] who found the WMAP results to 
be consistent with the concordance ACDM model. The statistic is limited by the one- 
sky realization we are able to measure, also called cosmic variance. To circumvent this 
problem it is suggested to utilize the fact that the CMB polarization is sourced by 
the local temperature quadrupole. On the one hand, the polarization signal contains 
information about the quadrupole moment at the reionization epoch. This could lead 
to a probability for the observed quadrupole moment of the order of 10~^ fH] compared 
with the concordance model. On the other hand, future measurements of the linear 
polarization of the CMB towards clusters of galaxies can betray the local quadrupole 
moment at the location of the cluster [201 1^ I22| thus circumventing the cosmic variance 
limit. 

Furthermore, we would like to mention that "the question whether the geometry 
of the three-dimensional space of astronomy might be non-Euclidean" has already 
been posed by Schwarzschild EHl llEj in 1900, fifteen years before the founding 
of general relativity! He stated the problem as follows (see p. 32 in (23)- "^s must 
be known to you, during this century [meaning the 19th century] one has developed 
non- Euclidean geometry (besides Euclidean geometry), the chief examples of which are 
the so-called spherical and pseudo-spherical spaces. We can wonder how the world 
would appear in a spherical or pseudo-spherical geometry with possibly a finite radius 
of curvature . . . One would then find oneself, if one will, in a geometrical fairyland; 
and one does not know whether the beauty of this fairyland may in fact be realized 
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in nature." Schwarzschild "actually estimated limits to the radius of curvature of the 
three-dimensional space with the astronomical data available at his time and concluded 
that if the space is hyperbolic its radius of curvature cannot be less than 64 light years 
and that if the space is spherical its radius of curvature must at least be 1600 light years" 
[SB]- In an Addendum to his article "On the permissible scale of curvature of space" 
P3 ] he mentioned also the possibility of spaces with non-trivial topology by referring to 
Clifford-Klein space forms. And he emphasized that such spaces do not necessarily lead 
to infinite universes as commonly assumed, even in the case of Euclidean or hyperbolic 
geometry. Schwarzschild concluded that the only condition imposed by astronomical 
data is that the volume of the Universe must be larger than the visible system of stars. 

It is the purpose of this paper to investigate two models of the Universe whose 
global topology is not the universal covering space of their spatial geometry. The spatial 
geometry of the models is the pseudo-spherical or hyperbolic space, i. e. it has negative 
curvature. The corresponding hyperbolic universes are non-compact and possess an 
infinitely long horn. The first model was introduced in 1976 by Sokolov and Starobinskii 
P7] and is described by a fundamental cell JF having infinite volume. The second model, 
called the Picard model jSHj, has also an infinitely long horn, but nevertheless possesses 
a finite volume. The Picard model has been investigated in detail in the context of 
quantum chaos [23 1^01 IHII • 

The Sokolov-Starobinskii model has been studied already in and |32] (see also 
the review In it has been claimed that the periodic horn topology produces a 
flat spot in the temperature map of the CMB even when multiple images of astronomical 
sources are unobservable. The fiat spot, i.e. a large flat region in the CMB sky map 
corresponding to a negligibly small metric perturbation, is argued in ^22^ to result from 
the exponential decay of the eigenmodes in the direction of the horn. 

We shall show in the following that the fiat spot found in [32] is a consequence 
of a too low wavenumber cut-off kc. For the cut-off chosen in |321, the eigenmodes 
cannot produce indeed a perturbation in the horn at the position which corresponds 
to the distance of the surface of last scattering (SLS). As a consequence, the sky maps 
computed in [32] do not show temperature fluctuations in the horn. However, if the 
cut-off is sufficiently increased, we shall demonstrate that there is no suppression in the 
horn and therefore no flat spot. 

The universes described by the Sokolov-Starobinskii and Picard model, respectively, 
possess negative spatial curvature in contrast to the concordance model corresponding 
to a spatially flat universe. Therefore the question arises whether a negatively curved 
universe is realized in nature. In |33t IM] we have analysed the CMB data and the 
magnitude redshift relation of supernovae type la in the framework of quintessence 
models and have shown that these data are consistent with a nearly flat hyperbolic 
geometry of our Universe if the optical depth r to the SLS is not too big. The restriction 
comes from the large amplitude of the fluctuations at large scales in the CMB. However, 
it was shown in [HSI that by replacing the quintessence component having a rest frame 
sound velocity of = c by a generalized dark matter component [IHl with a vanishing 
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rest frame sound velocity = 0, the amplitude in the CMB anisotropy is much smaller 
at large scales. Thus with such a generalized dark matter component universes with 
negative curvature are permissible even for larger optical depths r. 

Recently, an elhpticity analysis of the CMB maps has been reported |SZl EHl 1^1] 
which gives further support to a hyperbolic spatial geometry of the Universe. In these 
analyses hot and cold anisotropy spots in the CMB maps have been studied in terms 
of shape for various temperature thresholds. Analysing with the same algorithm the 
COBE-DMR, BOOMERanG 150 GHz and WMAP maps, an ellipticity of the anisotropy 
spots has been found of the same average value (around 2) for these experiments. The 
WMAP data confirm the effect for scales both smaller and larger than the horizon at 
the SLS. This suggests that the effect is not due to physical effects at the SLS, and can 
arise after, while the photons are moving freely in the Universe. 

Finally, we would like to mention that recently a model was presented t4Qj possessing 
positive curvature and a finite volume with the shape of the Poincare dodecahedral 
space. The authors of ref. PU| calculated the CMB multipoles for / = 2, 3 and 4, set 
the overall normalization factor to match the WMAP data at Z = 4 and examined the 
prediction for / = 2 and 3. They found a weak suppression of the power at / = 3 and 
strong suppression at / = 2 in agreement with the WMAP observations. However, in 
ref. jlO] only the modes up to kc = 30 have been used, and it thus remains the question 
about how this low cut-off affects the integrated Sachs- Wolfe (ISW) contribution. Our 
experience shows that increasing the cut-off usually enhances the ISW contribution. 



2. The CMB anisotropy 

The CMB anisotropy is computed for universes whose background models are described 
by the Friedmann equation (c = 1) 



a' \ ^ 8nG K 



" -Kl?) = — ^"-^5 • ^^{-1,0.+1} . (3) 

Here a = a{r]) denotes the cosmic scale factor as a function of conformal time V = f ^ 
and K the spatial curvature. The total energy density etot is given by a sum of 4 
components: radiation, baryonic matter, cold dark matter and a positive cosmological 
constant A. In this paper, we consider only scalar perturbations, i. e. vector and tensor 
modes are neglected. Furthermore, we are mainly interested in fluctuations on scales 
larger than the size of the horizon at the time of recombination since these fluctuations 
should contain the fingerprint of a non-trivial topology. Above the size of the horizon 
the temperature anisotropy can be computed neglecting physical processes important 
only on small scales such as the Silk damping. The metric with scalar perturbations is 
written in the conformal-Newtonian gauge in terms of scalar functions $ and as 

ds"^ = a'^{r]) {{l + 2^)dr]'^ - {1 -2'^)-fijdx'dx^} , 

where $ = \E' for an energy-momentum tensor T^^, with Tij = for i 7^ j and 
i,j = 1,2,3. {'jij denotes the metric of the 3-space.) The metric perturbation $ is 
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expanded into the eigenfunctions ipki^) of the (negative) Laplace-Beltrami operator A, 
i. e. (A + E)tpk{x ) = with k = y/E + K, e. g. for a discrete spectrum 

oo 

$(r/,f) = J2^kM^kA^) ■ (4) 

n=l 

In the single fluid approximation, the radiation, the baryonic matter and the cold 
dark matter are treated as one fluid which is valid on scales larger than the horizon size 
at recombination. There the entropy perturbation 6S is thus negligible, i. e. 6S = 0. 
Then the evolution of the metric perturbation $(77, x) gives in first-order perturbation 
theory in the conformal-Newtonian gauge 

^" + 3H{l + c1)^' - clA^ + {2H' + {l + 3cl){H^ ~ K)}^ = , (5) 

where H := a' /a. The quantity cj. = (S + l^mat/f^r) ^ can be interpreted as the sound 
velocity. Specifying the initial perturbation $ at = such that it corresponds to a 
scale-invariant Harrison- Zel'dovich spectrum (see jillj 17 = 2^/^TJl^^/^lmsi.t) 

$.(0) = ^ " and ^'.(O) = -M^ , (6) 

where $^(77) denote the expansion coefficients of $(77, x) (see equations (fT^ and (jHj) 
below), allows the computation of the time-evolution of the metric perturbation $. 
Here we have defined the density parameters fix '■= £x{vo)/£cTit with ^crit = 3Hq / {8ttG) 
for X = T (=radiation), mat (= baryonic + cold dark matter (CDM)), and curv (= 
curvature) with flcmv = 1 — ^tot- The Hubble constant Hq = h x 100 km s~^Mpc~^ 
is set in the following to h = 0.65. {a is a normalization constant and rjo denotes the 
conformal time at the present epoch.) The metric perturbation $ in turn gives the input 
to the Sachs- Wolfe formula which reads for isentropic initial conditions 



^(n) = 2$(r/sLS,x(r/sLs))- ^$(0,f(0)) + 2 T rf^^^MMl , (7) 

from which one obtains the desired temperature fluctuations 6T of the CMB. Here 
n denotes the unit vector in the direction of x(r7sLs) ~ ^ivo), i-e. in the direction from 
which the photons arrive. The single fluid approximation has been used in a lot of earlier 
papers concerning the CMB anisotropy for non-trivial topologies. We summarized this 
approach here, because we need it for the comparison with earlier papers. 

We use in our computations the tight- coupling approximation in which only the 
radiation and the baryonic matter is approximated as a single fluid before recombination, 
whereas the cold dark matter component is independent. The following summary of 
the tight-coupling approximation follows the lines of Mukhanov jlHI- The Sachs- Wolfe 
formula reads for the tight-coupling approximation 



^(n)= I d'k 



'?=^SLS 



+ 2 I d'k I dr]^^i;j:iTiri),9,^) 



dr] 
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where the fc-integration has to be interpreted as a summation in the case of a discrete 
spectrum. {T{7]),6,(j)) denote the spherical coordinates of the photon path in the 
direction n. Here, 5^^{ri) is the expansion coefficient of the relative perturbation in 
the radiation component, and V^^{ri) is the expansion coefficient of the spatial covariant 
divergence of the velocity field of the tightly coupled radiation-baryon components. 
The term in equation (jHJ which has to be evaluated at r/ = r/sLS is computed using 
the tight- coupling approximation. The second term, i. e. the integrated Sachs- Wolfe 
effect, is computed with a decoupled baryon component. In the case of isentropic 
initial conditions, the relation between the perturbation in the radiation and baryonic 
component during the tight-coupling phase is = f^j,- "^^^ time-evolution of the 
perturbations is determined by the following system of differential equations which has 
to be solved numerically. For the metric perturbation, one obtains {E^ = k'^ — K) 



-^^dv) + -Et + -T^- - K \ ^Ari) = 

dfj ''^'^ a'{r])\3 ^ \a{r]) J \ k\i) 3 / 



where the energy perturbation 6e^{ri) before recombination is given by 

(^) 

= ecDM(r/)%cDM(^) + , 
with c1 = (3 + jEh/sr:) ^, and after recombination by 

S^kiv) = ecDM(r?)%cDM(^) + eb(r/)%b(??) + e,{ri)5^/ri) . 
The energy and momentum conservation equations for the CDM component read 

^ (%CDm(^) - S-^-fcl^)) = -«(^)^fc,CDM(^) 

and ^ 

^ («'(^)^a;,cdm(^)) = aiv)E-,^^ir]) , 

respectively. For the radiation component, the energy conservation equation is 



whereas the momentum conservation equation is before recombination given by 



and after recombination by 
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After recombination the energy and momentum conservation equations for the baryonic 
component are analogous to that of the CDM component. The initial conditions of the 
various components depend all on that of '^'^(0) as 

= ^%b(0) = picBuiO) = -2<|.,K0) 

and 

%(o) = ^^,.(0) = i^,,cdm(o) = ■ 

The initial power spectrum $g(0) is taken as the usual Harrison-Zel'dovich spectrum 
having a spectral index = 1 which is in agreement with the current observations 
^1133]. Furthermore, we do not take a running spectral index into account since, as 
shown in only the first three multipoles C2, C3 and C4 are responsible for a non- 
vanishing running spectral index, and it is just these multipoles which are most strongly 
modified by a fundamental cell of finite volume. 



3. The Sokolov-Starobinskii model 

The model was introduced by Sokolov and Starobinskii |2Zj for cosmological studies. In 
order to describe the model, we have to introduce the hyperbolic three-space for which 
a wealth of models exists, see e.g. Here we choose the upper half-space model 

-^3 = {{xi,X2,X3) eR^ \xs>0} , (9) 

equipped with the Riemannian metric 

ds"" = ^^1 + + dxl ^^^^ 

X3 

corresponding to a constant Gaussian curvature K = —1. The geodesies of a particle 
moving freely in the upper half-space are straight lines and semicircles perpendicular to 
the xi-X2-plane. 

In the Sokolov-Starobinskii model one identifies points periodically along xi and 
X2 according to {xi,X2,X3) = (xi + /ia, X2 + z^&, X3) with /i, G Z. The positive constants 
a and b define the toroidal topology. The Sokolov-Starobinskii model has an infinitely 
long horn in the positive X3 direction whose physical cross section decreases according to 
(fTUj) with increasing X3. On the other hand, towards X3 — > 0"*" the cross section increases 
without bound. Thus the fundamental cell J-' = {x eM.^\0 < xi < a,0 < X2 < b,X3 > 0} 
has infinite volume. 

The eigenmodes of the Laplace-Beltrami operator on Sj^ 

A - x^ (— — —\ - X — 

consist of decaying modes and plane waves. The decaying modes are given by {k > 
0,(m,n) G (No x No) - {0,0}) 

, f sin^xi \ / sin^X2 A 

^kn,n{x) = NkmnXsKikiQXs) 2nm ^ nr.. ^ ' ^^^^ 

\ COS ——Xi I \ COS —r-X2 I 
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with 



N, 



kmn 



2k sinh nk (2 — 6mo) (2 — (5„o) 



ab 



and Q = 27r 



72 



> 



Here K[k{z) denotes the modified Bessel function of order ik. Due to the above 
normahzation factor N^rnm the eigenmodes (lllj) are normahzed as follows 



dlJ^i'kmnix)tpk'm'n'ix) = 5{k - k')5m,m.'5n,n' with dfl 



dxidx2dx3 



Xn 



Three examples of the eigenmodes (jllj) for k = 20 and three different pairs of m and 
n are shown in figure El The domain, where the modes dominate, decreases in xs with 
increasing Q. 

The metric perturbation corresponding to these modes can then be expanded as 



$(r/,f ; 



dk ^kiv) X] ^(^)"^)^) '4'kmnix] 
m,n 



(12) 



where the prime at the sum indicates that the term with m = = is excluded. The 
r{k,m,n) are real Gaussian expansion coefficients with mean value zero and variance 
independent of k such that the amplitude ^k{v) describes the usual power spectrum 
This leads to the usual statistical properties of the fiuctuations far away from the horn. 
The coefficients r{k, m, n) do not carry any information about the direction of the horn. 
The information about the horn is only contained in the eigenmodes of the Laplacian. 




Figure 2. Three eigenmodes using two sines for the eigenvalue fc = 20 and 
a — b — 1 are shown. The eigenmodes belong to {m,n) = (1,1), {m,n) = (2,1) 
and (m, n) = (3,4), respectively, from left to right. With increasing values of Q, 
the exponential suppression of the if-Bessel function occurs earlier in the horn. The 
suppression towards X3 is due to the a;3-factor in eq. (|11|) . 



The eigenmodes flllll decay exponentially in the horn due to the K-Bessel function. 
In addition to (fTT|) there are further eigenmodes which constitute non-decaying plane 
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wave solutions in the horn. They are given by {k > 0) 

,s,c.^^ 1 / sin(fclnx3) \ 

(x) = ^ x-s , X • (13) 

^ V cos(A;lnx3) ^ 

These plane wave solutions are normalized similar to ()11|) as 

j f^/^C(^)^fc'(^) = 5{k-k') . 

The metric perturbation can be expanded analogously to p2|l . but without the m and 
summation 

POO 

$^•^(77, x) = / cifcr(fc) $fc(r/)^^'^(x) , (14) 
Jo 

where r{k) are real Gaussian expansion coefficients. 

These solutions were not taken into account in the computation of the density 
perturbations in j23 1^2] • In [IZj it was argued that they should be omitted since 
they are not bounded as X3 — > cx). Constructing small perturbations from (fT!?jl would 
then require a fine-tuned superposition of these solutions, which would contradict the 
statistical hypothesis. 

Here, a remark is in order concerning the numerical evaluation of the fc-integrals ()12p 
and ()14|1 in the case that the integrands contain a Gaussian random function r{k,m,n) 
or r{k). In the case of a discrete /c-spectrum there is no problem. The integral is then 
replaced by a sum, as in equation (jlj), and each mode /c„ occurring in the spectrum 
comes with an independent Gaussian random variable. But in the case of a continuous 
/c-spectrum there arises the question of how small the step size has to be in order to 
approximate the integral numerically. Since the integrand is erratic, one cannot rely 
on the usual arguments concerning an approximation as, say, a Riemannian sum. Even 
infinitesimally close values of k possess uncorrelated coefficients due to the Dirac-delta 
6{k — k') in the correlation function {r{k)r{k')) . E.g. consider a unit interval Ak = 1 
and approximate the integral on Ak by terms. It is clear that in the limit N ^ 00 
the approximation to the integral approaches zero due to the increasing number of 
cancellations among the random terms. This implies that the normalization factor a 
which is used to normalize the fluctuation to the amplitude of COBE or WMAP depends 
on the discretization used in the approximation of the integrals. The value of has 
at least to be so large that there are many evaluation points on an interval where the 
integrand without the Gaussian random variable is nearly constant. Consider the mode 
ipl{x), eq. (fT^. having the simple Xa-dependence X3 sin(/c Inxs). This implies that in 
order to simulate a Gaussian random field there is for a fixed a maximal value of X3 
such that there are enough evaluation points for each sine-oscillation. Thus one has to 
ensure that there are no parts of the SLS so high in the horn, i. e. have corresponding 
high values of X3, that the CMB simulation is not valid for the chosen value of A^. We 
use in the following simulations a Gaussian quadrature with 16 evaluation points on the 
unit interval Ak = 1. This is sufficient for the chosen cosmological parameters. Finally, 
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let us remark that, if there are several contributing integrals, their normalization factors 
can be different if their integrands without the random variable oscillate on different 
spatial scales. 

As a first application let us consider only the decaying modes (fTT|) as it is done 
in [2Z1IS2]- In the temperature fluctuations 6T/T are computed for the Sokolov- 
Starobinskii model and it is claimed that the periodic horn topology produces a fiat 
spot in the sky map of the cosmic microwave radiation. The fiat spot, i. e. a negligible 
perturbation in the horn, is claimed to be the result of the asymptotic behaviour of 
the eigenmodes (jlll) which are exponentially declining for z = Qx^ > k, because of 
Kik{z) ^ \/^^~^ i'^ this range. Modes with k < 211x3 give a negligible contribution 
at 0:3 where the smallest possible value for Q, i.e.Q = 2tt in the case a = 6 = 1, is 
considered, since these are the modes which reach as far as possible into the horn. In 
[32 the cut-off in k is so low that the modes cannot produce a perturbation in the horn 
at the position which corresponds to the distance of the surface of last scattering. As 
a consequence, the sky maps computed in do not show temperature fluctuations in 
the horn. Increasing, however, the cut-off from the low value k^ = 10 used in |32.^ to 
the value kc = 140 shows that there are indeed fluctuations in the horn as it will be 
demonstrated below. 




50 100 150 200 250 



Figure 3. The fc-dependence of ^k{ilShs)/^k{0) in the single-fluid approximation is 
shown for a model including radiation with fimat — 0.3 and Q\ — (full curve) and a 
nearly flat one with fimat = 0.3 and Ha — 0.65 (dashed curve). 

There is nevertheless a subtle effect which can still suppress the amplitude of the 
perturbations in the horn even in the case of a high cut-off kc- As the above discussion 
has emphasized, it is important to include modes with sufficiently high values of k. If, 
however, the time evolution of ^k{v) leads to a too strong suppression of these modes at 
the surface of last scattering, i. e. if ^^(^ls) declines as a function of k too fast at the 
values of k which generate anisotropy in the horn, then these modes cannot contribute 
as much as they would for a nearly fc-independent ^kivsLs)- A pure matter model is 
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considered in [32] for which the time evolution of $fc(?7)/$fc(0) is independent of k 

5(sinh^ 1] — 3r] sinh r/ + 4 cosh?] — 4) 



^mat 
k 



iv) 



^mat 



(0) 



(15) 



(cosh rj — 1)^ 

Therefore this effect does not occur in the calculation in [52]. In general, however, 
especially if one takes also radiation into account, equation (0) leads to a /c-dependent 
time evolution of $fc(?7)/$fc(0) as can be seen in figure 01 where $fc(?7sLs)/$fc(0) is shown 
for two models containing also radiation, one with f2mat = 0.3 and = 0, and a nearly 
fiat one with flmat = 0.3 and Q\ = 0.65. The model with Q^ot = 0.3 has a much stronger 
/c-dependence as the nearly fiat model. Notice that for a pure matter model one has 
^k{VsLs)/^k{0) = const. The value of ^k{v) Vsls plays an important role because it 
describes the metric perturbation at recombination and via eq. the CMB anisotropy. 




Figure 4. The CMB anisotropy ST for the Sokolov-Starobinskii model for fimat — 0.3 
and flA = with a cut-off kc = 10 using a Gaussian superposition of the modes Hll|l. 
The white "curves" belong to ST = whereas increasing darkness corresponds to 
increasing deviations from ST — 0. 



Let us now discuss the cosmic microwave anisotropy of the Sokolov-Starobinskii 
model, which is here computed using the tight-coupling approximation. In contrast to 
|32] , we are not using a pure matter model and include the naive and integrated Sachs- 
Wolfe effect as well as the Doppler effect. The observer is located at {xi,X2,X3) = 
(0,0,1). For the coefficients r{k,m,n) Gaussian random variables are chosen. To 
illustrate the dependence on the numerically chosen cut-off in k, consider at first a 
model with fimat = ^tot = 0.3 as in jS^]. In figure El the CMB anisotropy is computed 
using the very low cut-off at kc = 10. The coordinate system is chosen such that the 
horn lies at the equator at one fourth from the left. One observes a suppression of the 
anisotropy in the horn as discussed above. In figureElthe cut-off is increased to kc = 140 
and no suppression is observed. The figure El represents the anisotropy 6T along the 
equator in the coordinate system of figures |31 and El (the horizontal line bisecting the 
sky map). The horn lies at ^ = |. The full curve shows the anisotropy for a cut-off 
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Figure 5. The CMB anisotropy 6T for the Sokolov-Starobinskii model for flmat — 0.3 
and JIa = with a cut-off kc = 140 using a Gaussian superposition of tlie modes Hll|l . 
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Figure 6. The CMB anisotropy 5T along the equator of figures 01 and [S] The 

suppression in the direction of the horn at 6 — ^ is clearly visible in the case of the 
kc — 10 cut-off (full curve). In the case of the kc = 140 cut-off (dashed curve) there 
are now fluctuations towards the horn. 

at fcc = 10 revealing clearly a suppression of the anisotropy around 6 = ^. The dashed 
curve presenting the kc = 140 cut-off displays anisotropics also in the horn. 

Up to now we have computed the CMB anisotropy solely with the eigenmodes (|TT|) 
and have ignored the plane wave solutions p3|l . These solutions contribute significantly 
towards the horn. As can be observed from figure El there are fluctuations for a 
sufficiently high chosen cut-off in k, but these fluctuations are smaller than farther 
away from the horn. This is due to the current restriction to the eigenmodes (fTTj) . 
The plane wave solutions have significant contributions especially in the direction of 
the horn as shown in figure [7| We now compute the CMB anisotropy for the model 
with f2tot = 0.3 and a cut-off k^ = 140 using all three types of modes. We use for the 
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Figure 7. The CMB anisotropy 5T along the equator of the sine (full curve) and 
cosine (dashed curve) plane wave solutions H13|l using a Gaussian superposition. The 
contribution from the cosine plane waves is multiplied by a factor of — ^ in order to 
get fluctuations comparable to the sine plane wave contribution. 




Figure 8. The CMB anisotropy 5T for the Sokolov-Starobinskii model for fimat = 0.3 
and r^A — with a cut-off kc = 140. A Gaussian superposition of all three types of 
modes is shown, i.e. the modes and the sine and cosine plane waves H13|) . The 
contribution from the cosine plane waves is weighted by a factor of — |. 



coefficients r{k,m,n) Gaussian random variables. The result is presented in figure |H1 
showing the CMB sky map and in figure IHl showing the anisotropy along the equator. 
As can be seen in figures El and El taking the sine and cosine plane waves into account 
increases the amplitude of the fiuctuations towards the horn. Thus the fluctuations do 
not betray the horned topology. Since we use Gaussian random variables as coefficients, 
there is also no fine tuning in contrast to what is stated in fI7\ . 

After having shown that the horn can be masked by the CMB fluctuations, let 
us now turn to the statistical properties of the CMB anisotropy. Expanding the 
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Figure 9. The CMB anisotropy is shown at the equator of the sky map from figure 

IHI 

temperature fluctuations with respect to spherical harmonics Y/^ yields the expansion 
coefficients aim which in turn lead to the multipoles Ci, see eq. ((H). The angular power 
spectrum 6T^ = /(Z + l)Q/27r is shown in figure ITUlfor the model where the superposition 
of all modes is used, i. e. for the anisotropy shown in figure |H1 Only large angular scales, 
i. e. multipoles with / < 65, are shown since higher values of / would require modes 
above our numerical cut-off at = 140. The values of ST^ display a fiat spectrum 
up to / ~ 40 as it is also revealed by the WMAP data ^21- (The error bars of the 
WMAP data shown in figure ^1 do not include the cosmic variance since we compare 
one-sky realizations.) The large fiuctuations of 6T^ in our model arise from the fact 
that we compute these values for a fixed observer and not, like CMBFast or CAMB, as 
a statistical average. The WMAP data display a very low quadrupole moment which 
is not reproduced by this simulation. This is due to the infinite volume of the Sokolov- 
Starobinskii fundamental cell. This contrasts to the Picard cell which is discussed in 
the next section. 

Let us now turn to the correlation function C which emphasizes the large angular 
scales and thus the low /-range. In figure ITT] Cft^) is shown (full curve) in comparison 
with the WMAP data (dashed curve). A large deviation from the WMAP data is 
revealed near = 0° since our computation takes only the values of Ci with I < 65 into 
account. The model does not describe the correlation hole for -i? > 160° as observed by 
WMAP which is due to the very small observed quadrupole moment C2. 

Up to now we have only discussed a model with Qtot = 0.3. However, current 
cosmological observations now point to a nearly fiat universe. Therefore we show in 
figure [121 angular power spectra ST^ for models with f2tot = 0.9 and f^tot = 0.95. 
We use here only the eigenmodes (|TT|l up to kc = 140 and ignore the plane waves. 
In order to obtain the correct values for higher values of / the contribution of modes 
above k = 140 is taken into account by assuming that their spherical expansion with 
respect to the observer point yields Gaussian expansion coefficients aim, i-e. that they 
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Figure 10. The angular power spectrum ST^ — l{l + \)Ci/2tt is shown for the sky 
map from figure |H1 as dots. The cruxes with error bars are the first year WMAP data 

m 




Figure 11. The temperature correlation function C^ff) is shown as a full curve for 
the sky map from figure |H1 The corresponding WMAP curve is displayed as a dashed 
curve. The dotted curve represents the best ACDM model described in |I7] . 

are statistically isotropic. Under this assumption, the /c-integral in (jH)) is for k > 140 
evaluated analogously to the CMBFast method. In figure [T2K) we choose Qmat = 0.3, 
l^A = 0.6, in figureHSb) ^mat = 0.3, = 0.65 and in figure IT^) fimat = 0.35, = 0.6. 
In all three cases an approximately fiat spectrum is observed having up to values of / = 65 
fiuctuations of the same order as the WMAP data. The topological structure of the 
Sokolov-Starobinskii model produces no signature in the angular power spectrum. Thus 
if our Universe would possess a small negative curvature with the Sokolov-Starobinskii 
topology, the CMB anisotropy would not reveal this "horned" topology. 
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Figure 12. The angular power spectrum 6T^ = 1(1 + 1)Ci/2tt is shown (full dots) 
for three nearly flat Sokolov-Starobinskii models having a) fimat = 0.3, Ha = 0.6, b) 
^^mat = 0.3, r^A = 0.65 and c) flmat = 0.35, flA — 0.6. The cruxes with error bars are 
the first year WMAP data U^]. 



4. The Picard model 

The crucial difference between the Sokolov-Starobinskii model and the Picard model |2H] 
is that the volume of the Picard model is finite. This is achieved by cutting off the lower 
(xs —>■ 0) part of the Sokolov-Starobinskii model which leads to the infinite volume. In 
the Picard model one chooses a = 1 and 6=1 for the periodical identification in the Xi- 
and X2-directions, respectively. In addition to these two identifications, in the Picard 
model all points are identified with those which are inverted at the unit sphere having 
radius 1 and its center at a; = (0, 0, 0). To describe this transformation let us represent 
the points x G i^s by Hamilton quaternions, q = Xi + 1x2 + j^s, with the multiplication 
defined by i^ = — l,j^ = — l,ij -|- ji = plus the property that i and j commute with 
every real number. The inverse of a quaternion q = qi + qi^ + ijsj + q^^ 7^ is then given 
by q"^ = \q\~'^{,qi — Q'2i — ^aj — Q'4ij), where |gp = ?i + Q'2 + Q'i + Qa- The Picard group F 
is then generated by two translations and one inversion, 

q^q + 1, gH->g + i, q^—q^^ . (16) 

The transformations 7 in ^3 are given by linear fractional transformations 

q ■J q = {aq + b){cq + d)~^ ; a,b,c,dEC ; ad — be = 1 . 

The group of these transformations is up to a common sign isomorphic to the group 
of matrices of PSL(2,C). The inversion q —q~^ bounds the fundamental domain 
by the unit sphere from below. From the three generating transformations (fT^ one 
can compose the transformation q 1— >■ igi which corresponds to a rotation by 180° 
around the x^-axis. This rotation together with the inversion q —q~^ are the new 
identifications in comparison with the Sokolov-Starobinskii model. The fundamental 
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Figure 13. The fundamental domain of the Picard group and the identifications of 
the surfaces. 



domain of standard shape is then given by [2HI HE] 

= |g = xi + ix2 + ja;3 - ^ < xi < ^, < X2 < ^, \q\ > l| , (17) 

see also figure ^1 The fundamental domain is a hyperbolic pyramid with one vertex at 
oo and the other four vertices in the points Pi = + P2 = | + -^3 = | + 
and Pi = —\ + \i + ^j. The volume of this non-compact orbifold r\i^3 is finite 



where 



vol(r\i53) = = 0.30532186... (18) 



4 

i/eZ[i]-{o} 

is the Dedekind zeta function, and Z[i] = Z + iZ are the Gaussian integers. 

Due to the additional inversion, the non-trivial eigenf unctions of the Laplace- 
Beltrami operator are more complicated than in the Sokolov-Starobinskii case and are 
not known analytically. The eigenfunctions are the so-called MaaB waveforms jSU] which 
are automorphic, i.e. satisfy i'i'jq) = i'iq) V 7 G F, g G and therefore periodic in 
Xi and X2- It follows then that they can be expanded into a Fourier series with respect 
to xi and X2 

M^) = u{x,)+ Yl «wa:3i^ifc(Qx3)e2-'('"^i+"^^) , (19) 

(m+in)eZ[i]-{0} 



where Q = l-n ymF^rn^ and 

u{x3) = <^ ^ ^ ^ . (20) 

I ^23^3 + ^z^z In X3 if = 

(Here we write '^{x^ = '?/'(xi, X2, 0:3) instead of 'ipiq)-) Kik{z) denotes again the K- 
Bessel function whose order is connected with the eigenvalue E hj E = k"^ + 1. If a 
Maafi waveform vanishes in the cusp, 

lim ^r^ix) = , (21) 

X3— »oo 
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it is called a MaaB cusp form. MaaB cusp forms have b^^^ = b™^^ = b^^ = b'^^ = 
and are square integrable over the fundamental domain, {ip'^;!^^^, 'ip')^^^) < oo, where 

(V',V'')= / dix^P^' with dfi = '^''^'^Y''' ■ (22) 

The coefficients a^^^ in (|Tn|l are not free expansion coefficients as in the Sokolov- 
Starobinskii model, but are uniquely determined for each eigenvalue k. This is due 
to the additional symmetry operation of the Picard cell. 

According to the Roelcke-Selberg spectral resolution of the Laplacian [HIIE^IESI, its 
spectrum contains both a discrete and a continuous part. The discrete part is spanned 
by the constant eigenfunction ip'f^^^ = vol(jF)~^/^ and a countable number of MaaB cusp 
forms ipk^'^^y "^k^^^y '^k^^^j ■ ■ ■ which we take to be ordered with increasing eigenvalues, = 
Eo < El < E2 < E3 < . . .. The smallest nontrivial eigenvalue is Ei = 43.8522464. . . 
PU] . Recall that the corresponding modes (fTT|) in the Sokolov-Starobinskii model possess 
a continuous spectrum due to the missing inversion symmetry. 

The continuous part of the spectrum > 1 is spanned by the Eisenstein series 
ip^^^'^'^^x) which are known analytically [^ES]. Their Fourier coefficients are given by 

r Eiscn 

1 Ai^(l + iA;) 

tEisen 

1 Aa'(1 - ik) 

° 0F|AA(l + ifc)| ' ' y^\AKil + ik)\ ' ^ ^ 



and 6f '^^'^ = bf''''' = and 



Eisen 
^kmn 



^|A/c(l + i*)l 



E 

A,/^SZ[i],A/i=m+ira 



A 



(24) 



where 



Ax(s) = A7i-'Tis)CK{s) (25) 



has an analytic continuation into the complex plane except for a pole at s = 1. With 
these coefficients the Eisenstein series are real and normalized (V'f'**™; '^f^^'^^) = S{k — k'). 
In the case of the Sokolov-Starobinskii model one has akmn = 0, and 60 and bi were 
chosen to obtain the real plane wave solutions (fTSj) . 

Normalizing the MaaB cusp forms according to {'ip™^^ , = 1, we can expand 

any square integrable function (p G L'^{r\S)-^) in terms of MaaB waveforms, |15] . 



[X) = 

n>0 



POO 

E<C^'''^)C^'(^) + / dk{^f-,<p)^f-{x) . (26) 

Jo 



Since the discrete eigenvalues and their associated MaaB cusp forms are not known 
analytically, one has to compute them numerically. We use the MaaB cusp forms 
computed along the lines described in (See PTII3^ for further references concerning 
the computation of the MaaB wave forms. For earlier computations see jHTIj.) 

In figures IT^ and IT31 we present two CMB simulations having Qtot = 0.95, the first 
having fimat = 0.30 and the second fimat = 0.35. Here only the cusp forms up to kc = 140 
are used. The observer is located in the upper half-space at x = (0.2, 0.1, 1.6). As in the 
case of the Sokolov-Starobinskii model no suppression of anisotropy is observed towards 
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the horn which hes at the equator at one fourth from the left, as in the previous sky 
maps. The figures IT^ and IT7I show the corresponding angular power spectra ST^. These 
spectra have a very small quadrupole moment as it is observed by WMAP. The next few 
moments are also well suppressed. The fc-summation in (jH} is carried out for k < 140 
using the cusp forms. The contribution of modes with k > 140 is approximated assuming 
statistical isotropy as in figure IT^ using the density of modes as it is given by Weyl's law 
for the cusp forms [23 ■ The figures fTHl and IT^ demonstrate that both models describe 
the temperature correlation function C{'d) much better than the concordance model. 
Both models display a very small correlation at large scales, > 60°, as observed 




Figure 15. The CMB anisotropy ST for the Picard model for f^mat = 0.35 and 
r^A = 0.60 with a cut-off fee = 140 using the cusp forms only. The observer is located 
in the upper half-space at x — (0.2,0.1, 1.6). 
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Figure 16. The angular power spectrum ST^ = l{l + l)C;/27r is shown as dots for the 
sky map displayed in figure [T^ The cruxes with error bars are the first year WMAP 
data. 
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Figure 17. The same as in figure [TBI but for the sky map displayed in figure [131 



140 / 



by WMAP. In addition, the two models agree even at smaller scales, ~ 10°, with 
the observations much better than the concordance model. This is due to the finite 
volume of the Picard cell as can be seen in comparison with figure ^2 for the Sokolov- 
Starobinskii cell which has infinite volume. Figure ITT] shows that the theoretical curve 
for the Sokolov-Starobinskii model oscillates around the curve of the concordance model. 

In figure 1201 the CMB anisotropy 6T is shown for an observer who is much 
higher in the horn at a; = (0.2,0.1,5.0). This is indeed very "high" in the horn as 
can be seen by comparing the volume V] of the fundamental cell (|17|) "above" the 
observer in the direction of the horn, i.e. the volume with > xf^^, with the volume 
Vi = vol(r\i33) — V| "below". Due to the hyperbolic volume element dfi, the volume 
above the observer is V| = ^ J^bs dx^/x^ = \{xf'^)~'^, xf"^ > 1. In the previous 
simulations for the Picard model, the observer was located at ^3*^^ = 1.6 for which 
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Figure 18. The temperature correlation function C{i3) is sliown as a full curve for the 
sky map displayed in figure [T^ i. e. with ilmat = 0.3 and fl\ — 0.65. The corresponding 
WMAP curve is shown as a dashed curve. The dotted curve represents the best ACDM 
model described in |T7] . 
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Figure 19. The temperature correlation function C(i?) is shown as a full curve for the 
sky map displayed in figure [T51 i. e. with fimat = 0.35 and JIa = 0.6. The corresponding 
WMAP curve is shown as a dashed curve. The dotted curve represents the best ACDM 
model described in [T7] . 



one obtains Vj ~ 0.207 and V| ~ 0.0976, i. e. roughly one third of the total volume lies 
in the direction of the horn. For the observer of figure EDI high in the horn with x'^^^ = 5 
one obtains Vi ~ 0.295 and V| = 0.01. Thus this observer point is very unprobable but 
we want to stress here that even such an extreme position does not betray the horn 
topology in the CMB anisotropy. The figure 1^ displays the corresponding temperature 
correlation function C{^). Again there is very low correlation for > 20°. In this case 
the observed anticorrelation for -(9 > 160° is not well reproduced. But nevertheless the 
overall agreement with WMAP data is much better than for the concordance model. 
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Figure 20. The CMB anisotropy ST for the Picard model for fimat = 0.3 and 
r^A = 0.65 with a cut-off kc = 140 using the cusp forms. The observer is located in the 
upper half-space at x = (0.2, 0.1, 5.0). 
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Figure 21. The temperature correlation function C{'d) is shown as a full curve for 
the sky map displayed in figure |^ i. e. with the observer sitting "high" in the horn 
and for fimat = 0.30 and J^a = 0.65. The corresponding WMAP curve is shown as a 
dashed curve. The dotted curve represents the best ACDM model described in |47| . 

Now let us consider a superposition of the contributions of the cusp forms and the 
Eisenstein series. The following model has fimat = 0.3 and = 0.65. In figure 123 
the CMB anisotropy 6T is shown where the Eisenstein series is weighted by a factor 
a = I relative to the contribution of the cusp forms. This factor depends on the chosen 
discretization of the ^-integration as discussed above. Using a discretization finer than 
= 16, i. e. 16 evaluation points per unit interval, would allow a larger factor. Since 
the summation of the cusp forms is not arbitrary, there are values of for which this 
factor can be one. The corresponding temperature correlation function C{i)) shown in 
figure ESI displays a very small correlation for 'd > 60° and the observed anticorrelation 
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Figure 22. The CMB anisotropy ST for the Picard model for flmat = 0.3 and 
r^A = 0.65 with a cut-off kc = 140 for a superposition of cusp forms and Eisenstein 
series (a = The observer is located in the upper half-space at a? = (0.2, 0.1, 1.6). 
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Figure 23. The temperature correlation function C(i?) is shown as a full curve for 
the sky map displayed in figure 1221 using cusp forms as well as the Eisenstein series 
(f^mat = 0.30 and Q\ — 0.65). The corresponding WMAP curve is shown as a dashed 
curve. The dotted curve represents the best ACDM model described in |47|. 

for -(9 > 160° is matched very well. 
5. Summary 

A large part of this paper was devoted to the question of the existence of flat spots 
in the CMB sky maps for universes with a horned topology. By the example of two 
such universes, the Sokolov-Starobinskii model and the Picard model, we showed that 
the infinitely long horn does not lead to fiat spots, i. e. to a suppression of the CMB 
fluctuations in the horn, if the wavenumber cut-off kc is chosen sufficiently large. The 
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flat spots reported earlier j32l E] are thus seen as the result of taking not enough modes 
into account in the expansion of the metric perturbation. We conclude that the CMB 
sky maps do not reveal a signature of the horned topology in form of flat spots. 

Another main point of this paper was to show that a universe with a horned 
topology, but with a finite volume, such as the Picard model, can explain the loss 
of power at large angular scales in the CMB anisotropy, as observed by CORE and 
WMAP. However, if the volume is infinite, as in the Sokolov-Starobinskii model, we have 
demonstrated that the low quadrupole moment is not reproduced. There is, however, 
some controversy of how serious one has to take this low value. The small values of 
the first few multipole moments Ci have, on the one hand, not been taken so serious 
to require some new physics, but rather have been considered as a manifestation of 
cosmic variance or an unsufficient consideration of Galactic emission ^Hj or caused by 
the local supercluster [nZj. From this point of view the Sokolov-Starobinskii model is 
a viable model. On the other hand, the suppression has been taken as a hint to new 
physics such as new information on the inflation potentials IHH t IHH]. However, these 
potentials have to be fine-tuned such that the arising power spectra are suppressed 
around the present day cosmological horizon [Snj- Such a fine-tuning is not necessary in 
the case of a non-trivial topology in a universe with negative spatial curvature due to 
the Mostow rigidity theorem [SIlin2]- Thus if the curvature scale is fixed by the densities 
fix, all side-lengths of the fundamental cell J-' are determined, and in turn the comoving 
wavenumber k at which the spectrum is suppressed, not because of the initial power 
spectrum but simply because of the absence of modes below the lowest wavenumber ki. 

In models with a non-trivial topology, certain points at the surface of last scattering 
can be identical due to the periodicity condition. These matching points are located on 
pairs of circles with the same radius. Along two such circles the temperature fluctuations 
6T produced by the naive Sachs- Wolfe effect are the same, which is called the circles-in- 
the-sky signature [HSI- In |E1] a search for such circles in the WMAP data was carried 
out for nearly back-to-back circles, i. e. for circles whose centers have a distance greater 
than 170° and whose radii are greater than 25° on the sky. They found no signature 
and rule out all topologies having such circles. However, it should be kept in mind that 
only the naive Sachs- Wolfe contribution leads to identical temperature fluctuations. The 
integrated Sachs- Wolfe contribution arises on the photon path to the observer, which 
is not identified for the observer and the "copy" of the observer. The same is valid for 
the Doppler contribution, since the observer and its copy see another projection of the 
velocity, in general. Furthermore, there are a lot of other secondary contributions to the 
temperature fluctuations. Nevertheless, it is claimed in jBl] that the naive Sachs- Wolfe 
contribution is strong enough for the identification of circles-in-the-sky. The Picard 
model is not ruled out by this work, since there are no nearly back-to-back circles. For 
the model with l^mat = 0.3 and = 0.65 with xs = 1.6, we obtain 40 pairs, where the 
largest distance of the centers is at 145° which is not covered by the study in [HI]. The 
model with fimat = 0.35 and Qa = 0.6 has only 32 circle pairs. The number of paired 
circles increases if the observer is posited higher up in the horn, i. e. at larger values of 
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X3. For the extreme position at X3 = 5 one obtains 275 circle pairs with distances up 
to 168.6° for fijnat = 0.3 and = 0.65. These circles could have been detected in 
However, for more generic observers, which are not sitting extremely high in the horn, 
the separation of the circles is too small such that the model cannot yet be ruled out. 

In conclusion, we would like to emphasize again that the Picard model studied 
in this paper is in nice agreement with the observed suppression of power on large 
scales in the angular power spectrum of the CMB (see figures ITHl fT^ 1^ and EH j) . This 
is in contrast to the concordance model which does not reproduce the experimentally 
observed suppression at {} > 60° and the observed correlation hole at ^9 > 160°. If 
future observations will confirm the WMAP data but with smaller errors, this can be 
interpreted as a clear hint to a non-trivial topology of our Universe having negative 
spatial curvature and a finite volume. 
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